Analytical solutions to single scattering of SH waves by a cylindrical 
fiber with radially gradient interphases 


Jun Zhang*®”, Lian Liu’, Longhai Zeng’, Jiading Bao® 


“College of Aerospace Engineering, Chongqing University, 
Chongqing, 400044, PR. China 
’Chongqing Key Laboratory of Heterogeneous Material Mechanics, Chongqing University, 
Chongqing, 400044, PR. China 
“School of Mechanical and Electrical Engineering, Guilin University of Electronic Technology, 
Guilin,541004, PR. China 


Abstract: In this work, analytical solutions to the single scattering of SH waves by a 
cylindrical fiber with two specific radially gradient interphase layers are supported 
based on the method proposed by P.A. Martin (Martin 2002, JASA). In the first case, 1) 
shear modulus u(r) = e?4" and the square of wave number, k? , is a linear function of 


I/r, and 2) u(r) = e~6r* and k? isa linear function of r2. For example, analytical 
solutions to single scattering of SH waves by a Sic fiber with the above two interphase 
layers embedded in aluminum are presented. The calculated scattering cross sections 
are compared with values obtained from an approximate method (dividing the 
continuous varying layer into multiple homogeneous sub-layers). The two methods 
yield similar results. The contribution of this work benefits the validation of various 
numerical methods used in the inhomogeneous media. 
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1. Introduction 


The multiple scattering of waves occurs throughout our daily lives — the scattering 
of sound by water drops in fog, for example, and the scattering of light by fine particles 
in the air. Research on the multiple scattering of waves has a lengthy history by virtue 
of its ubiquity !. The multiple scattering of elastic waves in composites is still an active 
research topic, as these materials are relatively new. Due to multiple scattering effect, 
even if the composite’s matrix and inclusions are both elastic, the elastic waves 
propagating in the material remain attenuated; this is called attenuation caused by 
scattering. To date, there have been many theoretical models proposed to evaluate this 
attenuation, as accurately evaluating it facilitates the dynamic characterization and non- 
destructive assessment ” of composite materials. 


The extant research can be roughly split into two groups. The first group consists 
of methods based on the wave function expansion and configurational averaging 
technique, which was pioneered by Foldy *’. In these methods, the exciting and 
scattered waves of each inclusion in a composite material are expressed as series of 
wave functions. The expansion coefficients of the scattered waves for each inclusion 
are related to the corresponding values of the exciting wave through a matrix T, the 
components of which are determined by the boundary conditions of the corresponding 
single scattering problem 8. For SH waves, T becomes a scalar since there is no mode 
conversion. After a sequence of mathematical operations, a homogeneous system for 
the ensemble averaged expansion coefficients of the exciting waves or scattered waves 
can be obtained. Solving this homogeneous system yields the attenuation coefficient. 
The second group consists of studies on the self-consistent method ™!?, where it is 
assumed that each inclusion behaves as an isolated one in a medium with the effective 
properties of the composite. The exciting wave acts on each inclusion is the coherent 
wave. The attenuation coefficient is then obtained through various consistency 
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conditions; for example, the mean wave field equals the coherent wave field !? or the 
total scattering by the inclusions embedded in the effective medium effectually vanishes 
11 


It is clear that, both the above introduced two groups of theoretical models 
necessitate solving a corresponding single scattering problem $, that is, the scattering of 
a plane incident wave by a single inclusion embedded in an infinite matrix. The single 
scattering problem is solved to obtain the far-field scattering amplitude, f (0), in 
models proposed by Waterman and Truell °! for example, and is solved to obtain the 
scattering cross section. In most of the existing theoretical models, the fibers/ 
particulates are usually assumed to be perfectly bonded to the matrix while the 
interphase layer between the fibers and matrix is neglected. The interphase layer should 
be accounted for, however, as it features gradient changed properties that are generated 
via chemical reaction and atom diffusion during the manufacturing process (or created 
artificially to improve compatibility between the fibers and the matrix). 


Up to date, the effect of the functional gradient interphase layer on the attenuation 
of elastic waves in composite materials has been sparsely considered in several works 
14-18. In these previous studies, in order to solve the corresponding single scattering 
problem, the inhomogeneous interphase layer is usually divided into multi-layers with 
homogeneous properties to approximate a continuous varying one. The wave fields in 
the fiber, intermediate sub-layers, and matrix are then expressed as series of wave 
functions. The continuity conditions of displacements and stresses at all interfaces are 
listed and the resulting coupled linear equations for the expansion coefficients are 
solved. When the intermediate sub-layers are sufficiently thin, this yields accurate 
solutions '°. This treatment of the gradient interphase layer is straightforward and can 
be used to deal with any type of radially gradient profile, but the coefficient matrix of 
the resulting linear system is sometimes ill-conditioned. Additionally, for composites 
with a high contrast of properties between fibers and the matrix, the interphase layer 
must be divided into a large number of sub-layers in order to obtain convergent results. 
This increases the computational cost as well as the condition number of the coefficient 
matrix |”, 


The transfer matrix method is another candidate for the above interphase problem 
'8 In this method, the displacements and stresses at the inner surface of each 
intermediate sub-layer are related to the corresponding values on the outer interface 
through a transfer matrix M, the components of which are functions of the material 
properties and geometry uniquely defined by the sub-layer. Using the transfer matrix 
for each intermediate sub-layer, the displacements and stresses on the outer surface of 
the whole interphase layer are related to the inner surface. Assuming the fiber, 
interphase layer, and matrix are perfectly bonded, a linear system for the expansion 
coefficients of waves in the fiber and matrix is established and solved. For SH waves, 
the transfer matrix M is a matrix of size 2 x 2. Therefore, the final resulting linear 
system to be solved is of size 2, which is much smaller than that in the method described 
above. Similarly, the final coefficient matrix of the resulting linear system can be ill- 
conditioned if the sub-layers are thin. 


Though several approximate methods to solve the single scattering problem with 
a radially gradient interphase layer have been proposed, to date, analytical solutions are 
still elusive. In this work, analytical solutions to single scattering with a radially 
gradient interphase layer of two specific profiles are reported for SH waves. The 
derivation process follows the work of P.A. Marin 7°, in which general solutions to the 
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single scattering of acoustic waves by an inhomogeneous sphere with spherically 
symmetrical properties was investigated. The two transformations used in the current 
derivation process differ from Martin’s due to the different governing equations for 
distinct waves. We also demonstrated how to extend the proposed transformations to 
other waves. 


The remainder of this paper is organized as follows: Section 2 derives the 
governing equations for SH waves in a radially gradient medium. General solutions for 
two specific radially gradient materials are then derived in Section 3. Section 4 presents 
analytical solutions for the single scattering problem of SH waves by a cylindrical fiber 
with the two specific interphase layers and a detailed example. Section 5 provides a 
brief summary and conclusion. 


2. Governing Equations 


For SH waves, only the displacement in the out-of-plane has value and can be 
expressed as w = w(x,t), where x = (x,y) is the position vector and £ represents 
time variable. The corresponding stress components can be expressed as follows: 


Ow Ow 
Oi, =F = Ox, = 01, = 0,0,, = u(x) ax 1053 = u(x) T (1) 


where p(x) is the shear modulus, which is a function of position. Substitution of the 
stress components to the governing equations of elastodynamics yields the following 


equation: 
O Ow) 0 Ow Ow 
uZ) uZ) © 
where p(x) is the mass density. After a simple mathematical operation, Eq. (2) can be 
re-written as follows: 
1 1 ôw 
u(x) è (x) a 
where V and V?are the 2D gradient and Laplace operator, respectively. c(x) = 
4 UCx)/pCx) is the wave speed, which is also a function of position. 


Define w(x, t) = uô (x)u(x, t), where A is a constant which is determined to 
force u(x,t) into the following form: 


Vu(x):-Vwt+Vow= (3) 


1 O7u 
V-u+Ku = TAD 
c Ot 


(4) 


For SH waves, A can be set to be —1/2 and then K = T Val SHV . The 


values of A should be changed for other scalar waves. In this work, only the radial 
inhomogeneity has been considered; that is, u(x)and p(x) are functions of only the 
radial coordinate r, so K can be simplified as follows: 


n" ' t 2 

E E ee (5) 
2u 2ru 4\u 

, hence u(x) satisfies the following 
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For time-harmonic problems, u(x,t) =u(x)e 
equation: 

Va+(k(r)+K |7=0 (6) 
where k(r) = w/c(r). Next, we seek the solution to Eq. (6) in the following form: 


u(x) =u,(r)Y, (0) (7) 
where n is an integer and Y„(0) is a cylindrical harmonic function, such as the 
Legendre function P,(cos6). 


We have the following relation: 
V? (uY, )=u, V" (¥,)+2(Vu, )-(VY, )+¥,V" (u, ) (8) 


because u, isa function ofr and Y, isa function of variable 0, so (Vu, ) . (VY, ) =0. 


We also know that r”Y„(0) is a separated solution to the Laplace equation in the polar 
coordinate system. To this effect: 


2 
VY, =-“{Y,, n=(-2,+00) (9) 


2 n°? 
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By substituting Eqs.(7) and (9) into Eq. (6), it is clear that u„(r) should satisfy the 
following equation: 


2 27] 
al! ities BEKO u,=0 (10) 
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which is a second-order differential equation for u,(r). 


3. Two Cases of Radially Gradient Interphase Layers 
Two kinds of radially gradient interphase layers are considered here. 


Case 1: In this case, the shear modulus is an exponential function ofr and k?(r) 


is a linear function of r~* as follows: 
. 20 
u) = me, [KOY = + —— (11) 
where 4, B, k?, and a are constants that are specified according to the following 
relation: 
Play 7 Hr Hy = Hm (12) 
2 2. 72 _ 42 
k - =k;; k l =k; 


where Uf, Um and kf, km are the shear modulus and wavenumber of the fibers and 
matrix, respectively; a and b are the inner and outer radius of the interphase layer. After 
a mathematical operation, the expression for K(r) in this case can be expressed as 
follows: 


Ky =-p°-£ (13) 
r 
Then the corresponding Eq.(10) becomes: 
2 1 p) _ 27 
AN E csc de (14) 
ar ž rdr r r | 
In order to solve this equation, we make the substitution: 
u(r) =rw(x) with x=6dr4 (15) 


where c and d are constants. Here, c is set to be -5 and d=1. The general process to 


specify these two constants is provided in the appendix. 6 is a parameter can be 
selected at the operator’s disposal. Eq.(14) can be transformed as follows: 


2 2 _ 2 
ki B és 12a-B n ul nogi 
ô ô x x 


(16) 
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According to the relative values of (k? — £7), solutions to Eq.(16) can be classified 
into three types. 


Case 1(a) (k? — p? < 0) 
Here, we choose (k? — B7)/5? = -* and set k = (2a — B)/6 so that Eq.(16) 


becomes: 


w"(x)+ [i+ ee G —n’)x* |w, (x) =0 (17) 


which is known as Whittaker’s Equation 7°. The general solution is: 
w, (x) =4,M,.,(x)+ BW...) (18) 


n 


where M,,(x) and W,,(x) are Whittaker functions. 


Case 1(b) (k? — p? = 0) 
We chose (2a — p) = 46 andset v=n -5 so that Eq.(16) becomes: 


w" (x) + ix —v(v+1)x” |w, (x) =0 (19) 


which is related to the Bessel’s equation. The general solution of Eq.(19) is: 
w, (x) = Vx ÍA, Jap, (V5) + B,Y,,(Vx)} (20) 
where J,() and Y,() are the Bessel function of the first and second type, respectively. 
Case 1(c) (k? — B? > 0) 
We chose (k? — B*)/6% = 1 and set n = —(2a—)/26 and v=n -7 so 
that Eq.(16) becomes: 


w"(x)+[1-2nx"'-v(v +x? |w,(x) = 0 (21) 
which is the Coulomb wave equation. Its general solution is: 
Wa (x) = AF an (n, x) + BG (n, x) (22) 


where F,(ņ,x) is the regular Coulomb wave function and G,(7,x) is the irregular 
Coulomb wave function. 


Case 2: In this case, the shear modulus is a Gaussian function of r and k?(r) is 
a linear function of r? as follows: 


Br? 2 

u(r) = me” , [kK] =k + yr (23) 
where u, B, k?, and y are constants that are determined in the same manner as in 
Case 1. The expression for K(r) is: 


K(r) =28 - p'r (24) 
So the corresponding Eq.(10) becomes: 
Tie lde lasap- Pe u,(r)=0 (25) 
d rdr r | 
We make the similar substitution as Case 1 to solve this problem: 
u,(r)=r w(x) with x=6r? (26) 


Similarly, ô is a parameter that can be selected at our disposal. Eq. (25) is transformed 
to: 


teal A 1 kK +26 n- 


n" + 
w(x) | 48 48 x 4x? 


| (x)=0 (27) 


As before, according to the relative values of (y — £7), solutions to Eq. (27) can be 
also classified into three types. 


Case 2(a) (y — B? < 0) 
We chose (y — B7)/52 = —1 and set k = (k? +28)/46 so that Eq. (27) 
becomes: 


2 
w" (x) -1+ Kx r-e w (x)=0 (28) 
the general solution to which is: 
W, (x) = AM nn (x) + BW ni (x) (29) 
Case 2(b) (y — B? = 0) 


We choose k? +28 = 6 and set v = =- so that Eq. (27) becomes Eq. (19), 
the general solution to which is: 


w, (x) = Vx f4, (Nx) + B,Y, (vx)} (30) 


Case 2(c) (y — B? > 0) 
2 _ 2 _ k2+2B CAN n-1 
We chose y — p^ = 46° and set n = eee and v= = 80 that Eq. (27) 
becomes Eq. (21). The general solution of Eq. (27) is thus: 
W, (x) = AF a-n (n, x) a B,Gnayn (n, x) (31) 


4. Analytical Solutions of Single Scattering Problems 


h interphase with / 
Ba radial gradient properties 


= 


Fig. 1 Schematic of single scattering problem with a radially gradient interphase layer 


In this work, the single scattering of SH waves by a cylindrical fiber with the two 
interphase layers discussed above was investigated. Figure 1 shows a schematic 
diagram of the single scattering by a cylindrical fiber with a radially gradient interphase 
layer. The inner radius of the interphase layer is a and the outer radius is b. The incident 
wave is a plane SH wave propagating in the positive x direction with a unit magnitude. 
Under the exciting of this incident wave, the total wave in the matrix, interphase layer, 
and the fibers can be expressed as follows: 


dei’, (K,,r)+4,H\ (k,r)]cos(n0), r>b 


n=0 


w(r, 0) = ¥169.(") r)+D,Q,(r) |cos(nd), a<r<b (32) 


YB, J,(k, r)cos(n@), r<a 


n=0 
where £ọ = 1 and ¢, = 2 (n #0), HOO is the nth order Hankel function of the 
first type; Q, (r) and Q, (r) are two functions with detailed expressions dependent 
on the type and material properties of the interphase layer. In Case 1, @Q, (r) and 


P, (r) can be expressed as follows 


M,,, (Sr) £iW,,(Sr) |, k -p <0 


p0) =—[ 
1,2 Jrut) K,n K,n 
pa(r)= qo” Sr H (LSP), i? — B? =0 (33) 


1 
P12 (r) = ag). 


In Case 2, Q, (r) and p (r) are: 


F (n,ôr)+iG ,(n,6r)|, Ke -B?>0 
aes ng 


1 
= — = or’ )+ ôr’ -B° <0 
P(r) Te eaa lr) EW, aalr], v-e < 


y—B’ =0 (34) 


pa) = 
eS ru) 
1 
=——— dr°)+iG or’) |, -B° >0 
P(r) Tp | Rath r)ti n101, > Y B > 


2 
Ayn, By, Cy, and Dp are the expansion coefficients to be determined by the continuity 
conditions of displacements and stresses across the interfaces at r = a and b: 


wW =a = r=a_ 
aw) _ aw (35) 
j OF |p at OF |a 


As the expansion coefficients An, Bn, Cn, and D, (n=1,2,---,n,,,,) are obtained, the 


far-field scattering magnitude ! and the scattering cross section * can be calculated as 
follows: 


f(O)= Şe, + cos (n0) 


n=0 


o, => 56,4 


k, n=0 
For the sake of illustration, the single scattering of SH waves by a Sic fiber 
embedded in an aluminum matrix with a radially gradient interphase layer was solved 


(36) 


nN 
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in this work. The material properties are listed in Table 1. Figure 2 plots the calculated 
o, changes with kma for the interphase types (Case 1 and Case 2). The results of the 
approximate approach (i.e., dividing the interphase layer into multiple homogeneous 
sub-layers) are also plotted for comparison. Both approaches yield similar results. 


Table 1: Material properties of Sic and Al 
Sic Pr (kg/m^3) (GPa) 


3181 188.1 
Al Pm (kg/m^3) Um(GPa) 
2706 26.7 


0.35 0.35 
-Æ Current theoretical solutions -& Current theoretical solutions 
0.3; |€ Multiple sub-layers 0.3} |€ Multiple sub-layers 
0.25 0.25+ 
Ee 02 g oz 
= = 
©” 0.15 ©” 0.155 
0.1 0.1 
0.05} | 0.05} 
Case 1 Case 2 
85s 075 1 125 15 175 2 8s o7 1 12 15 175 2 
k a k a 


Fig. 2 Comparison of o between current theoretical solutions and approximate 
method for interphase layer types (Case 1, Case 2) at several frequencies. Note: a=71lum 
and interphase layer thickness is A = 0.1 a 


5. Conclusions 


This paper reported the analytic solutions to the single scattering of SH waves by 
a cylindrical fiber with two specific radially gradient interphase layers, using the 
method proposed by P.A. Martin. In the first case, the shear modulus is an exponential 
function of rand k? isa linear function of 1/r; in the second case, the shear modulus 
is a Gaussian function ofr and k? is a linear function of r?. Besides SH waves, the 
extension of the current method to other waves has also been mentioned. 


An example of the single scattering of SH waves by a Sic fiber in an aluminum 
matrix was calculated and the scattering cross sections were compared against the 
values obtained by dividing the continuous varying interphase layer into multiple 
homogeneous sub-layers. The two approaches yielded similar results. The main 
contribution of this work is that it provides several benchmark solutions for 
inhomogeneous problems. 
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Appendix: General approach to make the second transformation 
Define u,(r) =r‘w,(x) with x = ôr’ , c and d are two constants to be determined. 
The first and second derivatives of u,(r) with regard to r can then be expressed as 
ul (r) =cr w, (x) +r" 'w (x) dd 
u"(r) =ce(c—1)r °w, (x) +8d (2c +d —1)r Ww (x) 407d? w" (x) 
Substitution of the expressions of u/(r) and u’(r) into Eq.(14), we can obtain 


cr °w (x) +d (2c +d)r wh (x) + 67d r w(x) +... 
3 A2 
G pin Jr ng(a)=0 — 
r r 


Therefore, if Eq.(A.2) is required to have the form of the confluent hypergeometric 
equation, the coefficients c and d should satisfy 

2c+d=0 (A.3) 
The term of the first derivative of w' (x) disappears. Then Eq.(A.2) becomes 


(A.1) 


2 2 

Edr w" (x)+ G pe e £ Je (x) =0 (A.4) 
r r 

Then d can be set to be 1 and c= -2 . Then Eq.(A.4) is simplified as 


wE E rae 2 Nn )=0 (A.5) 


